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Abstract 

A spin decoherence mechanism is proposed for localized electrons. The ir- 
regular phonon phase disturbances originated from phonon relaxation can 
influence electron spin precession with a net effect of spin phase decay. A 
quantitative analysis demonstrates relatively high efficiency of this mecha- 
nism in the low temperature and low magnetic field regime compared to the 
spin-flip processes. 
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I. INTRODUCTION 



Recently, much attention has been devoted to electron spin relaxation in quantum dots 
(QDs) since they provide a natural candidate for the qubits in quantum computing. A 
typical approach to this problem is to calculate the spin transition probability associated 
with the spin-flip processes, i.e., longitudinal spin relaxation. However, quantum computing 
is qualitatively limited by the processes that result in the destruction of electron spin phase 
coherence. For example, phase diffusion of localized electron spin can be characterized by 
relaxation mechanisms that are not related to spin-flip processes under certain conditions. 
Hence, further investigation of transversal (or phase) relaxation T 2 is crucial for accurate 
understanding. 

One such mechanism was proposed in Ref. 1 where a random change of spin precession 
and subsequent spin phase diffusion is associated with the transitions between electronic 
quantum states with different g factors. Although generally efficient, this process is frozen 
out at low temperatures due to its phonon-mediated nature and the direct spin-flip is ex- 
pected to be the dominant mechanism of phase relaxation. However, the spin-flip relaxation 
reveals a very strong (4th to 5th power) dependence on the magnetic field, 2 ' 3 becoming 
rather ineffective at low fields. Hence, it is necessary to explore other potential sources of 
decoherence, particularly in the low field and low temperature regime. In this work, we show 
that the spin-phonon interaction, which heretofore was considered mainly with respect to 
the resonant processes, can provide such a mechanism if a finite phonon damping is taken 
into account. 

Our analysis is based on the representation of spin-phonon interaction in terms of fluc- 
tuating effective magnetic field Q (in units of energy) acting on the electron spin "s > . This 
field is assumed to be composed of additive contributions Q p from each phonon p = {if, x} 
with a wave vector if and polarization x, i.e., Q = £ p f2 p . For the moment, let us focus 
on a single phonon contribution. Then, in the frame of reference rotating with the Zeeman 
frequency, the electron spin performs precession around the small Q p , which oscillates with 
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a phonon frequency u p . No alteration in the electron spin phase occurs due to such a har- 
monic perturbation with a possible exception of spin phase shift A0 O acquired at the initial 
period of interaction < t < lix joj v due to a random phonon phase 9 P . 4 

A different situation can be realized when a phonon harmonic oscillation is interrupted 
and resumes at a series of instant times t Vl and t 2 i (i — 0, 1, ...), respectively. The reason of 
such phonon fluctuations can be lattice anharmonicity, phonon scattering at the impurities or 
lattice defects, etc. These irregular phonon perturbations affect the electron spin precession 
resulting in the phase shift A0j at each interval of time tu — tu. Subsequently, the net effect 
of spin phase change <f> p (t) due to a phonon mode p can be expressed as <f> p (t) = SjA0j, 

(*2i < t). 

Note that for a large number of small changes A0j, their total effect can be described by 
a diffusion equation. Its solution leads to an exponential decay of electron spin phase with 
a relaxation rate T p x = \ (A</>?) r^ 1 , where r p is the mean time between sequential instants 
tu (or t 2 i) - 5 To estimate the spin phase change A0j caused by a phonon perturbation during 
the tu — tu, it is helpful to recognize that a single oscillator influence does not change a spin 
phase during its full period At p = 2n/u p as well as for any n integer periods n2n/u p . Hence, 
A0j can be approximated as a spin rotation Q p At p in an effective field Q p independently 
on duration t 2 i — tu- With the mean value (A0j) on the order of Q p /u p , one can expect 
Tp 1 ~ T^ilp/ujp for the phonon mode p and T 2 X ~ SpA^r" 1 ^/^ when the contributions 
of all phonons (with the population factor N p ) are taken into account. 

The qualitative consideration provided above shows that electron spin phase relaxation 
can be strongly affected by phonon phase damping of any origin such as phonon decay. 
Moreover, since this mechanism does not involve energy exchange, only the longitudinal 
(with respect to the external magnetic field ~B) component Q z of the effective fluctuating 
field is relevant to our case. These characteristics qualitatively distinguish the mechanism 
under consideration from other processes, most of which are determined by fluctuations of 
transversal components Q x and Q y at the resonant frequency with the Zeeman splitting. 
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II. THEORETICAL MODEL 



For a detailed quantitative analysis of the proposed mechanism, let us start with the 
spin-phonon interaction operator 

H s _ ph = Sis, (1) 

where the a-th component (a = x, y, z) of the fluctuating field takes a form linear in the 
creation and annihilation operators a) p and a p of the phonon mode p [—p = {— if, x}]; i.e., 

= £ w» = E v « W - a -v) ( 2 ) 

V V 

with a matrix element V£ of the spin-phonon interaction. The specific form of V£ will be 
discussed later. 

Now we focus on the spin evolution caused by random fluctuations of Q a . Obviously 
electron spin follows each of such fluctuations that result in its irregular behavior at the time 
scale r c of the Q a fluctuations. Actually a random single spin fluctuation associated with 
each phonon scattering is expected to be very small and drops out of the problem; instead, 
the total result of these small fluctuations averaged over the time scale At (r c <C At <^ T 2 ) 
is the subject of our investigation. The time evolution of mean spin value s can be described 
by the quantum kinetic equation 5 in the case of anisotropic medium and interaction H s - P h 
[Eq. (1)] 

±s(t)=ujxs(t)-rim-so], (3) 

where u is an effective field with components uj{ = J2j dij^sBj. 6 As usual, gij are the 
components of g tensor, the subscripts i and j relate to the crystalline coordinate system, 
Hb is the Bohr magneton, (...)= Tr{e~ Hd ^ T . . .} /Tre~ Hd / T where is the Hamiltonian of 
the dissipative subsystem (lattice vibrations in our case), and T is the temperature, uj and 
T are expressed in units of energy. The matrix T of relaxation coefficients is composed of 
Fourier transformed correlation functions 
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l, v = l,u H = (n„ (r) n„> w = ^J{% (r) n„> e^dr (4) 

— oo 

with f2 M (r) = exp(iH d r)Q^ exp(—iH d r) . It has a canonical form in the frame of references 
x, y, z with z directed along ul (so that /i, v — x, y, z). With a provision that the 
correlation functions are symmetrical, 7^ (uj) = 7^ (a;), the matrix T has a simpler form: 
F xx = 7T (7° z + n7 ro ), T ra = 7T (7° 2 + 717^), T 2Z = irn(-f xx + 7 W ), r M „ = -irwy^ , (^f), 
where 7 ° 2 = 7jm (0), n = n(w) = (l + e-/ T ) /2, s = -±ztanh(u;/2T), u = u z = ^^f" ■ 

One can see that the coefficients r xa; and T ra responsible for transversal relaxation consist 
of two parts, = ^n^ yy (or iin^ xx ) and T 2 ~q = 7i^ z . Comparison with the longitudinal 
relaxation coefficient T zz shows that the term T 2 ~^ stems from the contribution of spin-flip 
processes involving energy exchange between the Zeeman and phonon reservoirs. Since the 
longitudinal relaxation has been the subject of a number of recent studies, 2 ' 3 ' 7 ' 8 we focus on 
the analysis of T^q term. 

The correlation function Fourier image [Eq. (4)] of the effective field Q is expressed in 
terms of the phonon operators according to Eq. (2). In turn, the Fourier image of phonon 
correlation functions (p p (uj) = (Q p (t)Q~ p ) w is Lorentzian-like since the corresponding Green 
function satisfies the equation G p (uj) = (oUp/n) [uj 2 — uj 2 — 2uj p M p (uj)]~ 1 , where the "mass" 
operator M p {uj) depends on the phonon interaction (see, for example, Ref. 9). In the most 
general case, this correlation function takes the form 10 

= 1 (2Ne±11IM , (5) 

* - ^|)7^ 2 + r» 

where T p (uj) = ImM p (w) depends on the specific mechanism of phonon scattering. In such 
a manner, T p (uj) is a function of temperature due to the anharmonicity of the third and 
fourth order; furthermore, there are contributions by other sources of phonon scattering 
(point defects, isotopes, dislocations, crystal boundaries and interfaces) that reveal different 
dependencies on uj and u p . Hence, evaluation of the relaxation coefficients becomes too 
complicate to be approached analytically. Instead, to proceed further, we utilize the phonon 
relaxation time that can be extracted from the thermal conductivity measurements (see 



Ref. 11 and the references therein). An expression appropriate for the correlation function 
Fourier image 7^ (cu) was derived in Ref. 5 in the relaxation time approximation. For our 
particular case of u> = and /i — v — z, it can be reduced to 



where r p = 1/Y p {uj p ) is the relaxation time of phonon mode p (i.e., phonon lifetime). In 
most cases, one can assume oj p ^> r" 1 and neglect the second term in the denominator of 
Eq. (6). Then, along with the definition of the operator Q p [see Eq. (2)], one can express 
the non-resonant phonon contribution to the transversal spin relaxation rate in the form 



which is in accordance with the qualitative analysis discussed earlier in this paper. The 
phonon population factor N p is given as [exp (u p /T) — l] -1 . 

Equation (7) is the starting point of our investigation on the proposed spin relaxation 
mechanism. However, this still requires the detailed knowledge of the phonon dispersion uj p 
and the relaxation time r p for each phonon mode p. By taking into account the conditions 
frequently encountered in quantum computation utilizing semiconductor QDs, we restrict 
our consideration to the case when the radius a of the electron state is much larger than the 
lattice constant and the temperature is sufficiently low. Since the spin-phonon interaction 
matrix Kf is significant only for the phonon wave vector q < l/a , a large a essentially 
limits the summation of Eq. (7) to long wavelength phonons. Subsequently t p , which is a 
complex function of the temperature and phonon frequencies, 12 can be considered in the 
long wavelength limit. Moreover, at low enough temperatures T < T bs (T bs rs 10 K in the 
case of Ref. 12), only one term originating from the boundary scattering survives for phonon 
relaxation. 11 Since this mechanism is insensitive to the temperature as explained by Ref. 11, 
it is adequate to assume a constant phonon relaxation time r p ~ r p h for long wavelength 
phonons at T < T bs . This permits us to avoid the problems associated with the complex 
dependence of uj p and t p , which can be very specific for each particular sample. 
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A. Effect of g-factor fluctuation 



To evaluate VJ, we consider the spin-lattice interaction via phonon modulation of g 
factor. In general, the spin-lattice interaction Hamiltonian can be written in terms of the 
tensor A ijH : 13 > 14 > 3 



H s -ph — AjjkiUBBjSjUki, 

ijkl 



(8) 



where % is the strain tensor Uki averaged over the electron ground state \g) = ip g (1^): 
Uki = (g\iiki\g). By way of important example, we consider a z-directed magnetic field 
and a localized electron with the axial symmetry with respect to the z-axis. This reduces 
Eq. (8) to the form of Eq. (1) with Vt z = [(A 33 — A 3 i ) + A 3 i A] pB; here, A denotes 
the dilatation A = Ejtijj and the Voigt notation is adopted (A 33 = A zzzz , A 3 i = A zzxx , 
A 6e = A xyxy ). Then, the matrix element of the spin-phonon interaction takes the expression 



V? 



h 



1/2 



[(A 33 - A 31 ) e p z q z + 5„, L A 31 q] $ (t) VbB, 



(9) 



2pVu> p 

where p is the mass density of the crystal, V is the volume of the sample structure, l? p the 
polarization vector of the phonon mode p, x = L, T, and $ (if) = (g\ e l Q ' r \g). The spin- 
lattice relaxation rate in Eq. (7) can be calculated by treating the phonon modes based on 
the isotropic elastic continuum model with the longitudinal and transverse sound velocities 
cl and ct- Assuming the axial symmetry for the local electron center, i.e. $ (if) = $ (x, z) 
(x = qao/2, z = q z /q, the parameter ao represents the electron state radius as mentioned 
before), one can obtain 



^o 1 = r^m 



x- 



'ph 



coth 



r e ff \ c 3 / 

^x ) F T (x) + 4F L (x) coth ' ' ' 



(A 33 - A 31 ) 2 p%B" 
27r 2 hpc T al 

|V + c)V (x,z)dz, 

J z\l - z 2 )<5> 2 (x,z)dz, 



X 



dx, (10) 



F L (x) 
F T (x) 



(11) 



where r^ 1 is an average phonon relaxation rate, Tfff = hc^/ksao is the effective tempera- 
ture, and C = A 31 / (A 33 — A 31 ) = —1/3 if one assumes that the strain induced part of the 
effective g-tensor = Ti^^A^um is characterized by zero trace, i.e., A 33 + 2A 3 i = 0. When 
a is much larger than the lattice constant, the upper limit x max in the integral of Eq. (10) 
may be taken to infinity since $(~(f) restricts the actual phonon wave vectors to q < 1/ao 
as discussed above. 

Let us evaluate spin relaxation of a shallow donor with an effective Bohr radius Ob (=a>o) 
and $ (x, z) = (1 + x 2 ) 2 . Utilizing the constant phonon relaxation time approximation 
r p ~ T p h for T < Tf, s , the integral in Eq. (10) can be evaluated analytically 

where T u(±) = (16/15tt) T^ } . Note that Eq. (12) is obtained with if \\ [001]. In the case 
of cubic symmetry [where only two constants A m and A 33 = —2A 3 i in Eq. (8) describe 
the effect of spin-phonon coupling], an expression T 2 ~ 1 for an arbitrarily directed can 
be obtained in terms of the direction cosines / = B x /B, m = B y /B, n = B z /B. Our 
calculations show that this is achieved by multiplying the factor 

f0/B) = 1 + - l) P, (13) 

to Eq. (12); P = 3(/ 2 m 2 + m 2 n 2 +n 2 l 2 ) , < P < 1. One can see that the angular dependence 
of our mechanism does not result in zero relaxation under any direction of 1^ . Moreover, 
the directions along the principal axes ([001], etc.) can result in maximal relaxation, while 
the same directions sometimes forbid the spin-flip processes. 13 ' 15 

As an example, we consider a Phosphorus shallow donor in Si with a# = 1.8 nm. The 
phonon relaxation time can be extracted from the low temperature measurements of Si 
thermal resistivity 16 in terms of the theory developed in Refs. 17 and 11 {r ph = 2.4 x 
10~ 8 s). The spin-phonon coupling constants were estimated in the works of Refs. 13 and 
15. However, we believe that direct determination of coupling constants by means of EPR 
measurements of Si:P under an applied stress gives more reliable data. A corresponding 
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experiment was performed in Ref. 18, where the constant A 66 = 0.44 was found. Similarly, 
our estimation obtained A 33 = 0.31 and A 31 = —0.155 that gives T^ 1 = 1.3 x 10~ 4 s _1 at 
the magnetic field of 1 T and low temperatures T <C T„(±) — 10 K. 

In another important case of a Si shallow donor in Alo.4Gao.6As, the data on EPR under a 
uniaxial stress 19 provide rather strong spin-phonon constants of A 33 = 19.6 and A 31 = —9.8. 
This gives the estimation T^ 1 = 6.1 x 10 -2 s _1 and 6.1 x 10~ 4 s _1 for the magnetic fields 
of 1 T and 0.1 T, respectively, at T = 4 K under the assumption that phonon lifetimes are 
identical in these crystals. 

Similar calculations can be performed for an electron localized in a QD of L xy = 2ao 
in the lateral width and L w = eL xy in the thickness. Under the condition e < 0.1, an 
approximate formula takes the form 




where the fitting coefficients are b T — 1, &l = c|/c|, ct = 0.33 — 1.27e 2 , dr = 0.35 — 0.395e 2 , 
c L = 0.97 - 28.5e 2 , and d L = 0.40 - 3.76e 2 . 

Let us compare, as an example, spin phase relaxation caused by the phonon decay 
[Eq. (14)] with the spin-flip admixture mechanism (Ref. 2) in a GaAs QD with L w = 3 nm 
and L xy = 25 nm, assuming r ph = 2.4 x 10~ 8 s and A 33 = 19.6. For the relatively strong 
magnetic field of 1 T and T = 4 K, our mechanism and the spin-flip mechanism give 
^2,0 ~ 0-1 s_1 an d l^f 1 = ^2,w = 10 s_1 ' respectively, while for B — 0.1 T both mechanisms 
predict almost the same rate of ~ 10~ 3 s _1 . In lower magnetic fields, our mechanism prevails. 

B. Effect of hyperfine constant modulation 

The (/-factor modulation described in Eq. (8) is not the only possible mechanism of spin- 
phonon interaction. For an alternative process, let us consider the hyperfine interaction 
(HFI) of localized electrons with the nuclei: 

H hf = a h f^2\^('7^j)\ 2 ~fj't, (15) 
3 
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where is the HFI constant and / j is the nuclear spin situated at site j with the position 
Lattice vibrations near the nuclear equilibrium positions can lead to effective field 
fluctuations and, subsequently, the spin-phonon interaction. Taking into account the long 
wavelength phonons with respect to the mean internuclear distance ~ n i ^ 3 (rij is the nuclear 
spin concentration), the main part of this interaction for a typical nuclear spin configuration 
can be represented as in Eq. (1) with 

Tt = nJl(I+l)m/V QD a hf A. (16) 

Here, the unit vector n is directed along the effective nuclear field defined by Eq. (15) and 
Vqd = (/ |V ; (^)| 4 d 3 !*) 1 . Calculation of the phase relaxation rate for the case of a shallow 
donor results in the expression, which is similar to Eq. (12), 



WV, 1 /, . T 1 



+ ^ (IT) 



where the parameter 



I(I + l)nia 



a/ WhpV QD c\al (18) 
is independent on the magnetic field. In the case of an electron localized in a QD, one can 
find the approximate rate through an analogy with Eq. (14): 



T^o^^^clf + dl^, (19) 

where Chf = 3.7 — 68e 2 , dhf = 2.7 — 9.8e 2 , and e < 0.1. Numerical estimations provided for 
a donor in Si and GaAs in terms of Eq. (17) indicate inefficiency of this mechanism with a 
very long relaxation time (about 10 14 s and 10 8 s, respectively). Hence, this mechanism can 
be neglected in most cases. 



C. Two phonon process 



So far, we primarily considered the influence of phonon decay on spin phase relaxation via 
linear spin-phonon interaction as given in Eq. (8). Namely, the effect of phonon scattering 
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with an electron spin on phonon relaxation has not been considered (i.e., electron spin- 
induced phonon decay). The Hamiltonian of this process can be derived in terms of spin- 
two-phonon interaction Hf} ph = ^ Dijki mn ^BBiSjUkiu mn with the spin-phonon coupling 
constants Dij k i mn . Now the fluctuating effective field takes the form Q a = J2 PP > WP' p 'Q p Q p > 
(WP' P are the matrix elements of H K s _ ph ), so the correlation function Fourier image 7^ (u) 
[Eq. (4)] is expressed in terms of phonon correlation functions ((Q Pl Q P2 ) (t) Q P3 Q P4 ) w - Its 
calculation performed in a harmonic approximation leads to a simple expression 5(cu pi — 
uj p 2 )(^pi,p3^p4,p2 + $P2,P3$P4,pi)(2N pi N P2 + N P1 + N P2 ). Substituting this function for 7^ (u) 
and a parameter D for the dominant contribution among the coupling constants Dijkimn, 
the spin phase relaxation rate for the two-phonon process is given at low temperatures 
(T < hcT/kBdo) approximately as 

2 - a ~ 21p>4 \ hc T ) ■ [ U> 

Parameter D can be estimated as D = 3(g — 2)C 2 /Eg (g, C, and E g are the electron g 
factor, deformation potential and energy gap). 1 Numerical evaluation of Eq. (20) at low 
temperatures (T = 4 K) predicts a long relaxation time. In the case of GaAs at B — 1 T, 
one can find T 2j o ~ 3 x 10 5 s, which is too long to be of any experimental or practical interest. 



III. DISCUSSION 

To illustrate the significance of the mechanism under consideration, let us briefly survey 
the most important spin decoherence mechanisms reported in the literature: the HFI and 
spin-lattice interactions. In the presence of the HFI, an electron spin performs precession 
around the sum of the external magnetic field ~§ and the effective field ~§hf caused by the 
HFI. Dispersion of ~B*hf over an ensemble of QDs results in a relatively fast electron spin 
phase diffusion (see Refs. 20 and 21); however, it causes only a partial dephasing (< 67%) 
and can be essentially eliminated as B » B h f (B hf < 1 G for typical Si QDs). 

In the case of a single electron in a QD, the electron spin can change its phase through the 
HFI since the nuclei also perform precession around the effective field caused by the electron 
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spin. This field proportional to ^i/? j)\ 2 [see Eq. (15)] is inhomogeneous over the QD 
volume, which distorts the mutual correlation of nuclei spin configuration and subsequently 
causes an alteration in the direction and strength of i^/. 20 ' 22 However, this relaxation is 
rather long and can be suppressed if B » B h f. In addition, it can be further reduced in 
the case of full nuclei spin polarization 22 and/or isotope purification. Hence, the spin-lattice 
(i.e., phonon) interaction provides the most fundamental and unavoidable source of electron 
spin decoherence. 

Among the spin-lattice interaction mechanisms, the phonon-mediated transitions be- 
tween the ground and excited states modulate the precession velocity leading to very ef- 
fective decoherence, 1 when their energy separations are small enough. However, under the 
assumption ksT <C 5 this relaxation is reduced as exp(— <5o/A:bT). Thus, the spin-flip pro- 
cesses and the phonon-decay induced mechanism considered in this paper provide the main 
contributions at low temperatures. Moreover, these two mechanisms differ in the magnetic 
field dependence. When the magnetic field decreases, the spin-flip process yields to spin 
phase diffusion induced by phonon relaxation as mentioned above. The estimated mag- 
netic field strength for this cross-over (e.g., < 0.1 T) is well within the range of practical 
importance. 

IV. CONCLUSION 

We considered spin phase diffusion of a localized electron through anharmonic phonon 
disturbances. In contrast to the spin-flip process where only the resonant (with the Zeeman 
energy) phonons are relevant, electron spin phase acquires random shifts when relaxation of 
any (resonant or nonresonant) phonon occurs. A quantitative analysis shows that the con- 
sidered phase relaxation reveals a relatively weak dependence on the magnetic field strength 
and the temperature compared to the direct spin relaxation processes or other mechanisms 
that involve the excited electron states. In addition, a specific dependence on the magnetic 
field direction [Eq. (13)] is attributed to this mechanism. Thus, one can expect that at 
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low temperatures and magnetic fields the spin phase diffusion mediated by the phonon re- 
laxation can become dominant over the spin-flip processes. As for quantitative estimation 
of the relaxation rate, the decisive role belongs to the phonon lifetime T p h- In the present 
study, we estimated T p h from the experiments conducted in bulk Si. It is not apparent if 
this estimation is applicable to the case of QDs. Moreover, the phonon lifetime may be a 
function of geometry and composition of the structure under consideration. However, the 
qualitative signatures of the proposed mechanism is expected to persist and may provide a 
ground for experimental verification. It should also be pointed out that the framework of the 
developed theoretical model allows more accurate estimation when the detailed information 
on phonon dispersion and relaxation is taken into account. 
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